Prospects for transient gravitational waves at r-mode 
frequencies associated with pulsar glitches 
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Abstract. Glitches in pulsars are likely to trigger oscillation modes in the fluid interior of 
O neutron stars. We examined these oscillations specifically at r-mode frequencies. The excited 

^J ', r-modes will emit gravitational waves and can have long damping time scales C(minutes - days), 

j.^ ■ We use simple estimates of how much energy the glitch might put into the r-mode and assess 

&Q' the detectability of the emitted gravitational waves with future interferometers. 



O I 1. Introduction 

One possible source of burst-like gravitational waves is a pulsar glitch. These are anomalies in 
the observed spin frequency of a neutron star (NS) wherein there is a sudden increase in its 

^vq . on timescales of a few days-weeks. In most cases, these glitches are produced by young pulsars. 

It is believed that these glitches are a consequence of global events within the star's fluid interior 



* 



and its solid crust [T] . These glitches could excite a variety of quasi- normal modes [2] . A previous 
search for gravitational waves associated with a glitch in the Vela pulsar in 2005 focussed on 
/-modes whose frequencies are believed to lie in the range 1-3 kHz [3] and last O(100 ms). 

A possibility to be considered in gravitational wave searches of this type is the excitation 
of r-modes oscillations during a pulsar glitch [5]. These r-modes, which are a type of inertial 
modes (oscillation modes restored by the Coriolis force), have frequencies proportional to the 
star's angular velocity and they are considered relatively more efficient GW emitters compared 
to other inertial modes [5]. Also, these modes evolve in time with an exponential decrement 
that depends on the mode's frequency and dissipation processes [6]. Reference [7] details that 
when a rotational lag between the fluid component and the crust reach a critical value, an 
instability that excite r-mode oscillations within the NS interior sets in. Because of this lag, the 
fluid component (already 'pinned' to the crust by a series of vortices) will be unpinned and will 
transfer angular momentum to the crust which will consequently spin up. As a result of this 
process a pulsar glitch will be produced. In that work, the authors considered short-wavelengths 
modes at which the quadrupole moment I = m = 2 would not be excited. 

In section [2] we describe the form and duration of the GW signal one might expect from 
r-mode oscillations. In section [3l we explore the energetics of the r-mode signals associated 



with pulsar glitches and the plausibility of detecting GWs from r-modes. Finally, in section 2] 
we summarise our findings. 

2. Gravitational waves signals from r-modes in neutron stars 

We parameterise the expected GW signal from an r-mode oscillation in terms of harmonic 
ringdown gravitational waves signals [8], 

h+(t) = h cos(uj r t)e~ t/Tr h x (t) = h sm(uj r t)e~ t/Tr , (1) 

where h+ and h x are the gravitational wave signal polarisation 'plus' (+) and 'cross' (x), 
uj r = 2irf r is the frequency of the r-mode oscillation, r r is the damping time of the signal and 
ho is the amplitude of the signal at the detector given by equation (23) in [9J 

ho = J—~aoj*MR 3 J. (2) 

V 5 c° r 

This expresses the amplitude of the signal in terms of the NS mass M, radius R, distance from 
the observer r, a dimensionless parameter J (1.635x 10~ 2 ) that depends on the stellar mass 
distribution and the amplitude of the mode oscillation a. We consider the amplitude a in more 
detail in section [3l 

According to [ID], for sufficiently slowly rotating stars, the r-mode oscillation frequency f r 
is proportional to the spin frequency of the star / S pm- It is expected that the £ = m = 2 mode 
dominates the GW emission [11] . Under these assumptions the r-mode frequency is given by, 

4 

uj r « -Q, (3) 

where Q is the angular spin frequency of the star (= lixjspin) OH]. There are several sources 
of uncertainty in this relation. These include rotational effects [12] . which modify the mode 
frequency by an amount that depends upon the average density, and relativistic effects [IT] . 
which introduce fractional corrections of the order of the stellar compactness M/R, so that 
equation ([3]) is accurate only to the ~ 10% level. 
The damping time of the r-mode is given by, 
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- = + , (4) 

Tr ''"viscosity TQRR 

where T v i scos it y is the dissipation time-scale due to the viscosity in the fluid component of the NS 
and tgrr is the gravitational radiation reaction time-scale. Following (DJ, the radiation reaction 
time-scale is a function of the mode frequency, 

%tJMR\ (5) 



TGRR 15 2 C 7 

Note that the damping time of the signal r r is a balance between dissipative processes with 
time-scale r v i scos it y > and the radiation reaction time-scale tqrr < 0; while r v i scos i ty acts 
to dissipate energy and reduce the oscillation amplitude of the mode, tqrr forces the mode's 
amplitude to grow due to the emission of gravitational radiation [13]. In this work we assume 
that the mode is in the stable regime so that 7Viscosity<<| 7 "GRR|j and so the viscosity time-scale 
dominates in the duration of the signal. 

Levin and Ushomirsky proposed in [J2] the existence of a viscous boundary layer in the 
crust-core interface of the NS which yields the following estimate of the dissipative time-scale, 

1 =0 .0L,-^^y /2 (^Y, (6) 
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where -Rio, M1.4, /> and pb are the radius (10 Km), mass (1.4 M@), the density at the crust-core 
interface, and an estimate of this density (pb = 1.5 x 10 14 g cm -3 ), respectively. Tg is an assumed 
internal temperature of 10 8 K in the NS. The value of the fitting parameter F depends on the 
dominant scattering processes in the fluid interior [15] : 

1/15 for electron-proton scattering, 

F ~ < (p/pb) 5 ' 4 for neutron scattering, (7) 

5p/pb for electron scattering. 

Here, we assume that neutron-scattering dominates and we have set p = pb- 

Finally, the ratio 5u/u measures the fractional velocity mismatch between the crust and 
the core. As demonstrated in |14j . for a sufficiently slowly rotating star, the crust does not 
significantly participate in the r-mode oscillation, so that 5u/u ~ 1. Just how slow the rotation 
needs to be for this to be the case depends upon the thickness of the crust, but Figure 1 of |14] 
indicate that all young glitching pulsars are likely to have 5u/u ~ 1. Note, however, that in 
more rapidly rotating stars Su/u ~ 0.1, lengthening the viscous decay timescale by a factor of 
~ 100. 

3. Detectability 

The detectability of the GW signal from r-mode oscillations excited by a pulsar glitch is 
determined by the total amount of energy that the glitch deposits in the r-mode, and by the 
fraction of that energy which is radiated as GWs by the r-mode oscillations. 

Following [9], the instantaneous gravitational wave luminosity from a source at distance r, 
undergoing damped oscillations with initial amplitude ho, angular frequency u and damping 
time r is, 

£GW = ^rV(M^) 2 . (8) 

The total time-integrated energy emitted in GWs is, therefore, 

&E GW = C Q^ r2 ^ h l T - ( 9 ) 

The mode energy is 

E = a 2 n 2 MR 2 J. (10) 

A rough estimate of the total energy associated with a pulsar glitch is (see e.g. [3]), 

^glitch » m 2 (J^pj , (11) 

where I the moment of inertia of the NS, O is the spin frequency and AQ/O, is the size of the 
glitch relative to the spin-frequency. An upper limit on the mode amplitude a can be obtained 
by assuming all of the energy associated with the glitch is channelled into r-mode excitation 
and we find, 

a =b) (it) • (12) 

where / (typically ~ 0.261) is a dimensionless parameter dependant on the stellar mass 
distribution [6]. To evaluate the total power contained in short-duration, narrow-band signals, 
it is convenient to use the root-sum-squared amplitude: 
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h\{t) + h 2 x (t) dt 




h r 1/2 , (13) 



where the right hand side is the approximate value for the damped sinusoid we consider for our 
GW signal. Combining equations ([2J, (fTT|) and (fl2|). the root-sum-squared GW amplitude from 
r-mode excitation in terms of the stellar parameters and the size of the glitch is, 



2-kGM 3 f~ z AU > 1/2 
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Figure 1. Sensitivity curves of gravitational waves detectors and upper limit estimates h r 
the gravitational wave signal as a function of the r-mode frequency f r = uj r /(2n). 
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In figure [TJ we compare the estimated GW upper limits from glitch-induced r-mode 
excitations in various known pulsars with the noise spectral densities of the initial LIGO 
instrument during its fifth science run [16], the advanced LIGO detector [17j.|18j (currently 
under construction) and the proposed 3 rd generation Einstein Telescope (ET) [19] . We compare 
the estimates for frequently glitching pulsars (observed to have glitched more than five times, 
shown in crosses) and all other pulsars, shown in dots. The glitch observations are compiled 
from [2U1[2TJ[22] an d [23]. We find that, in this optimistic case, the estimated GW amplitude is 
a factor ~ 6 below the noise floor of even ET in its most favorable configuration 



4. Discussion 

We have presented some estimates of the energy associated with r-modes excited by pulsar 
glitches and the corresponding estimated upper limit GW amplitudes. We see that, even in 
the most optimistic scenario of energy transfer from the glitch to the mode, the GW amplitude 
is well below the sensitivity curves of both existing and planned gravitational wave detectors. 
Specifically, even for the most rapidly spinning pulsar we consider, where the r-mode frequency 
lies near the most sensitive part of the detectors, the expected amplitude is ~ 4 x 10 -26 . 

However, the strong frequency dependence of the r-mode oscillation implies that rapidly 
rotating NS have a significantly greater chance of detection in 3 generation instruments such 
as ET. In addition to this, it is possible that these faster rotators possess a significantly smaller 
slippage parameter 5u/u (see section [2]) , leading to longer duration and more detectable GW 



signals. There is clearly some interest in exploring r-mode excitation in more rapidly rotating 
stars. 

Related to this, we end by noting that, despite the greater frequency of glitches occurring 
in young pulsars, there is a report of a millisecond pulsar (MSP) glitch in |24j . As well as 
pushing the expected GW signal closer to the sensitive regions of interferometric GW detectors, 
higher spin frequencies yield intrinsically stronger signals. Accreting millisecond pulsars present 
a particularly tantalising opportunity for future analyses and the GW r-mode detectability from 
such sources is to be considered in future work. 
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